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1 Introduction 



We consider the two obstacles problem for monotone operators (possibly degenerate or 
singular) of the type 



with a Dirichlet boundary condition in a bounded domain C W 1 . 

The two obstacles problem for the operator (1.1) consists of finding u € such that 
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= {v € W 1,G (O) : V < v < ^ a.e. in Q, } 
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ip,<p £ W 1)G {VL) (for the definition of the Orlicz-Sobolev space W 1 ' G (fl) see next section), 
where G is related to a by 

G(x,t) = [ a(x,s)sds, x € Q, t>0. (1.4) 
Jo 

Wq' G (Q) is the closure of Cg°(n) in W 1 ^^) and is a Banach space of weakly differen- 
tiable functions v with 

/ G(x, \Vv(x)\)dx < oo. 

In particular, when in (1.1) we have a(x,t) = t p ^~ 2 , with p(x) > 1 a given bounded 
function in fi, we deal with problems involving variable growth conditions, the so called 
p(x)-Laplacians. The study of such problems has been stimulated by problems in elasticity 
(see [32]), in fluid dynamics (see [3], [9], [29], [33]), image processing models [8] and problems 
in the calculus of variations with p(x)-growth conditions (see [2], [22], [23], [32], [34]) and 
some more general class of differential operators (see [4], [6], [12], [33]). 

Here we are specially interested in the more general class of quasi-linear operators of 
Ladyzhenskaya-Uraltseva type (see [19], [20]), when a(x,t) : fl x M + — > R is given by a 
function measurable and bounded in x for all t > and Lipschitz continuous in t, a.e. x € fl, 
and, such that, there are positive constants a<a 

< a < fa fo^ + 1 < a for t > 0, (1.5) 
a(x,t) 

where at = da/dt, and also lim t _>.o+ ta(x, t) = 0. 
The assumption (1.5), in fact, implies 

(a(x, m ~ a(x, |C|)C) • (e - C) > 0, VC, C G M n , i + ( (1.6) 

and and lim^oo ta(x, t) = oo for a.e. x <G ft (see [7], for instance). 

As a consequence, we have the uniqueness of the solution and also the weak maximum 
principle for A. In this work, after recalling the natural functional framework of the Orlicz- 
Sobolev spaces associated with A, in section 2, we extend some classical properties on the 
bilateral problem (1.2) — (1.3) in this more general framework, for instance, by including 
examples like a(x,t) = a(x)t p ^~ 2 log(/3(x)t + 7(2;)) with bounded functions j(x), p(x) > 1, 
and a(x), f3(x) > a.e. in x € fi. 

In section 3, we use the continuity property of the truncation operator v t-^- v + = tiVO = 
sup(i>,0) for the strong topology of W 1,G (fl) to extend some continuous dependence results 
in W Q ' (Q) of the variational solutions to (1.2) — (1.3) with respect to the data. 

In section 4, we prove the Lewy-Stampacchia inequalities (here a A b = inf (a, b) and 
a V b = sup(a, b)) 

AcpAf<Au<A^jVf (1.7) 

in an abstract form, extending the approach of Mosco [25] to the two obstacles problem, that 
includes the above class of operators. Although the inequalities (1.7) are known, in particular, 
for linear operators (see [30]), our proof is new and more general. As a consequence, under 
additional (Holder) continuity hypothesis on 14 a(x, •), we obtain the same regularity for 
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the solution u of the two obstacles problem as in the equation without constraints (see [20] ) . 
For instance, for bounded obstacles <p, ip we conclude that u G C 1,a (£l), if we impose /, A<p, 
Ail) £ ^°°(^)) a regularity obtained in [19] with different assumptions. 

Finally, in section 5, we give two new applications to systems of obstacle type. In the 
case of the N- membranes problem, when u = (ui, . . . , un) has the constraint 

ui>U2>---> un a.e. in f2, 

we extend some of the results of [5], in particular, the C 1,a regularity and the strong approx- 
imation in (Wq ,G (Q)) N by solutions of a penalized system. For the case of a special class of 
implicit double obstacle problems, when the obstacles depend on the solution in the form 

\f(uj - tpij) <Ui< f\(uj + tpij), i = 1, . . . , N, 

where ipij, ipij are certain given positive constants, we are able to show the existence of a 
minimal and maximal solution for the corresponding system of quasi- variational inequalities, 
which is of the type arising in problems of stochastic impulse control (see [14], [25], [30] and 
[31]). 



2 Preliminaries on Orlicz-Sobolev Spaces 

The Orlicz spaces and the Orlicz-Sobolev spaces as well are defined for the Young functions 
(or N-functions) like G defined by (1.4). Let for all x £ Q, g(x, •) : R — > R is an odd, increasing 
homeomorphism from R onto R; g(x, t) > 0, when t > 0, while the function G : x R — > R, 



G(x, t) := / g(x, s)ds 
Jo 



for all x G 0, and all t > belongs to class G (see [26], p. 33), i.e. G satisfies the following 
conditions: 

(i) for all x G fl, G(x, •) : [0, oo) — > R is an increasing function, lim^oo g(x, t) = oo, 
G(x, 0) = and G(x, t) > whenever t > 0; 

(ii) for every t > 0, G(-,t) : SI — > R is a measurable function. 

For a Young function G, we define the (generalized) Orlicz class, 

Kg{&) = {u : O, —?■ R, measurable; pc{u) := l G(x, \u(x)\)dx < oo} 

Jn 

and also the (generalized) Orlicz space, 

L G (Q) = {u : Q -> R, measurable; lim / G(x,X\u(x)\)dx = 0}. 
which is a Banach space endowed with the Luxemburg norm 

\\u\\ G = inf j M > 0; J G^x,^l^x<l|. 
The (generalized) Orlicz-Sobolev space is defined as follows 

W 1 ' ^) = lu G L G {tt);^- G L G (f)),i = 1,2, ...,n|, 
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is also a Banach space with the norm: 

IMIi.g* = HIV^IHg + IMIg- 

These spaces are more general that the usual Lebesgue or Sobolev spaces, but many proper- 
ties of functions in these spaces can be extended. In particular, the Poincare type inequality 



/ G(x, \u\)dx < / G(x,c\Vu\)dx, 
Jn Jn 



holds for any u € W^' G ($7), where c is twice the diameter of (see [15]). 
The Holder inequality extends to (see [26], Theorem 13.13) 



/ uvdx 
Jn 



< C*|M|g|MIg> Vn G L G (Q),v e L G (n), 



where C is a positive constant, and G is the conjugate Young function of G, that is, 
G(x, t) = sup{ts - G(x, s); sG»}, ViGfi, t > 0. 

We also denote by G* the Sobolev conjugate of G, that is 

rt G _1 (x, s 



provided 



°°G-\x,s) 

as = oo. 



s (n+l)/n 

It is also well known, that 

L°°{n) ^ L G (n) ^ l 1 (o), 

with continuous imbeddings, and if £1 is a bounded domain with a smooth boundary, then 
the imbedding W ' (£1) L (Q.) is continuous. 

In this work given the function a : $7 x ]R + — > ]R + satisfying the assumption (1.5). Let 
the mapping g:$7x]R— ^Rbe defined by 

Then g satisfies the conditions (i) — (ii), and the corresponding G is a strictly convex Young 
function. 

In our case, (1.5) implies 

ta(x t) 

< 1 + a < ^7 ' ( < o + 1 ,a.e. x e tt t > 0, (2.2) 
G(x, t) 

and G satisfies the so called A2- condition (see, for instance [1]), which implies that L G {Q) = 
K G {n) (Theorem 8.13 in [26]), and that L G (Q) is a linear separable space. 

Relation (2.2) assures (see Proposition 2.2 in [23]) that L G (Vi) is an uniformly convex 
space and thus, a reflexive space. 
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In Orlicz-Sobolev spaces an important role is played by pc( u ) - the modular of the 
L G (£1) space. If u m , u <G L G (Q) then (see [23] for the proofs) when ||u||g > 1 it holds 
— Pg{ u ) < IKIIg" ) an( i we have 

\\u m \\G -> oo o pc(u m ) oo, 

and when ||«||g < 1 also ||u|| G +a < pc(u) < ||n||^ + -, which implies 

\\u m - u\\ G -> <S> /9 G (n m - u) -> 0. 

We refer to [10], [11], [18], [26] for further properties of (generalized) Lebesgue-Sobolev 
spaces. 

Here we need the following lemma: 
Lemma 2.1. If u G W^ G {VL), then u + ,u~ € W^ G {9.) and 

Du +={ Du i ^ u>0 ' and £>u" = ■[ °' ' /n -°' 
\ 0, z/u<0, \ ifu<0. 

Here « + = uV0, u _ = — uAO. This lemma holds in Wq ,G (Q) as well, and therefore these 
Orlicz-Sobolev spaces are closed with respect to 

u\Jv = u+(v — u) + and uf\v = u — (u — v) + . 

The proof of this lemma is due to Gossez (see [16]) and is basically the same as for usual 
Sobolev spaces (see, for instance, Theorem 1.56 in [30]). 

Lemma 2.2. The embeddings W^ +1 (n) ^ W 1,G (O) ^ W^~ +1 {n) are continuous. 

Proof. The second part was observed in [20]. In order to prove the first part, by Theorem 
8.12 in [1] it is enough to check, that t a+1 dominates G near infinity, which is true, since by 
taking the log in the right inequality of (2.2), we conclude that there exists c > and T > 
such that G(x,t) < ct u+1 for t > T.U 

Proposition 2.3. // u m -> u in W 1,G (Q), then u+ n ->• u + in W l ' G (£l). 
Proof. This follows by convexity of G and by the inequality (see [22] ) : 

G{x, t + s)< 2 1+a (2 + a) (G(x, t) + G{x, s)) a.e. a € fl Vt, s > 0. 

We have 

/ G(x, |V(u+ -u + )\)dx = j G(x,\x{ Um >o}Vu m -X{u>o}Vu\)dx 
Jn Jn 

~ C ( In X{Mm>0}G ^' ' V ^ m ~ U ^ dx + J n l x W>°} ~ X{ u >o}\G(x, \Vu\)dx) -> 0, 

and so Vu+ — > Vu + in L G (ft). Arguing in the same way, we get also u+ — > u + in L G (ft), 
which completes the proof. ■ 
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Remark 2.4. a) Assuming G(x,t) = G{t), i.e. G is independent of variable x, we say that 
L G and W l,G are Orlicz spaces, respectively Orlicz-Sobolev spaces (see [1]). 

b) Assuming G(x,t) = \t\ p ^ with p(x) G L°°(J7) ; p(x) > p > 1 a.e. in £1, we denote L G (Q) 
by L p ( x \VL) and W l ' G (VL) by W 1,p ^ x \Vt) and we refer them as variable exponents Lebesgue 
spaces, respectively variable exponents Sobolev spaces. 

c) Our framework enables us to work with spaces which are more general than those described 
in a) and b). Besides the example given in the introduction for a(x , t) = a 

7(x)) withp(x), j(x) > 1 and a(x), (3(x) > a.e. x G fi, we could consider any linear com- 
bination with positive coefficients or any composition of functions satisfying a condition like 
(1.5). 



3 Variational solutions 

1 G 

We introduce the energy functional J : W^{Q) -»• R by 



J(u) = [ G(x, \Vu\) dx, Vu G W^ G {Q) (3.1) 
Jn 



which is strictly convex, weakly lower semi-continuous and coercive in WQ' G (tt) (see [21]). 

1 G 

Moreover, J is Gateaux differentiable, and J'(u) = Au at u € W ' (Q) is given by (see [23], 
for instance) 

(Au,v) = [ a(x, |Vu|)V« • Vv dx, v G W 1,G (O). (3.2) 
Jn 

Hence Au G W' 1 ' ^) = (W 1,G (fi))', the topological dual of W 1,G (S1) and if we assume 

/ G L G (p.) C W~ 1,G (Q), G* being the conjugate Young function of the Sobolev conjugate 
of G, we can rewrite the problem (1.3) — (1.4) in the form: 
find u G K^, such that, 

(Au-L,v-u) >0Vt>GlK^, (3.3) 

where we set 

(L,v) = [ fv dx, Vv G Wq' G (CI). 
Jn 

Proposition 3.1. Under the condition (1.6) t/ie operator A is strictly T-monotone, i.e., for 
any u, v G W 1,G {£l) 

(Au-Av,(u-v) + ) > 0, (u- v) + G W 1,G (O). 

Proof. In fact, using Lemma 2.1, by (1.6) we have 

J | [a(x, \Vu\)Vu - a(x, \Vv\)Vv] ■ V(u - w) + jdx 
= / < [a(x, |Vu|)Vu - a(x, \Vv\)Vv] ■ (Vu - Vv) \dx > 0, 

J{u>v} I J 

if (u — v) + ^ 0, i.e. if meas{u > v } > 0.B 



In this section we assume (1.5) and 



y>, V G W 1,G (^) such that K% ^ 0, (3.4) 

1,G, 



for which it is sufficient to assume ip > ip a.e. in fi and both 99", V + G W„ 
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Theorem 3.2. The problem (3.3) has a unique solution u = u(f, (p, ip) G and is equiva- 
lent to minimize in the functional F, defined by 

F(v) = [ G(x, \Vv\) dx - [ fv dx, v G WQ' G (n). (3.5) 
Jn Jn 

Moreover, if u denotes the solution corresponding to f, <p, ip, then 

f > /, 9? > V?, ip > ip implies u > u a.e. in fi. 

Proof. The existence and uniqueness are standard results for strictly monotone, coercive 
and potential operators, as observed in more general Orlicz-Sobolev spaces already in [17]. 
The comparison property follows easily by the T-monotonicity (see, for instance, [25] or [27]): 
take v = u V u G K 1 ^ in (3.3) and v = u A u G in (3.3) for u. By addition, one finds 

(Aii - Au, (u - u) + ) + (L-L,(u- u) + ) < 0. 

Since L — L>0, and A is strictly T-monotone, one immediately deduces (u — u) + = 0, which 
means that u > u.M 

Remark 3.3. This argument also shows a weak maximum principle in the form: if Au > Au 
in SI and u > u on Oft in the sense (u — u) + G Wq' G (£1), then u > ii in ft. 

Similarly we have a "L°°-continuous dependence" property, even without the L°° regu- 
larity on the solutions: 

Proposition 3.4. 

\\u - n|| L oo ( Q) < \\ip - (p\\ L °°(ci) V - ^||loo(q), 

where u and u are the corresponding solutions of the problem (3.3) and (3.3) with the same 
f- 

Proof. Let I = \\ip — <p\\l°°(si) v IIV* — M\L°°(n) < °°- Set v = u + (u — u — l) + G in 
(3.3) and v = u — (u — u — l) + G in (3.3). By addition one gets 

h := / (a(x, \Vu\)Vu- a(x, |Vu|)Vu) -V(u-u-l) + < 0. 
Jn 

On the other hand, recalling (1.6) we know, that if meas{u > u + /} > 

h := / (a(x, \Vu\)Vu - a(x, |Vn|)Vn) • (Vn - Vn) + > 0. 

J{u>u+l} 

Since I\ = h, we conclude u — u — l < 0. Reversing the role of u with u we get u — u — l < 0.B 

Exactly as in Proposition 4.5 of [23], we have the following interesting result. 

1 G 

Lemma 3.5. Let (1.5) holds and u m — u weakly in W ' (Q). Then if 

limsup(Au m , u m — u) < (3.6) 

then u m — > u in Wg' G (f2) strongly, and Au m — > Au in W~ l ' G (Q) strongly. 
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Theorem 3.6. Under the assumption (1.5) let u rn be the solution to (3.3) with compatible 
data (f m ,(p m ,ip m ), such that 

fm—tf^ L G , ip m — > if and ip m — > i\> in W l,G (£l) strongly. 

1 G 

Then u m — > u strongly in W ' (Q), where u is the solution of the limit problem (3.3) with 
(/,</2,Y>)- 

Proof. For arbitrary v G we obtain that v m = (v A ip m ) V ip m G K^™ and, by 

1 G 

Proposition 2.3, v m — >■ v in VF" ' (0). Since, by coerciveness, ||« m ||i,G < C> there is u G 
W 1,G (f2), such that, for a subsequence u m — 1 u in W llG (ft) weakly and in L 1 (0). So n G 
and by lower semi-continuity we have 

F(u) < liminf F(u m ) < lim F(w m ) = F(v), Vt> G K?. 

m— >oo m— too ™ 

1 G 

Hence u solves (3.3) and then, since w m = (u A <^ m ) V tp m — s- u in Wq' (fi), 

limsup(An m , n m - u) < lim (Au m , u m - w m ) + lim (Au m , u? m - u) 

< lim / / m (w m - u m ) = 0. 

By Lemma 3.5, this implies the strong convergence of the whole sequence u m to uM 

Remark 3.7. In the case of only one obstacle problem, lower (ip) or upper (if) obstacle 
(corresponding to take formally tp = +oo or ip = — oo respectively) all previous results of this 
section hold in similar way. For instance, Proposition 3.4 takes the form 

\\u - «||L°°(n) < U - ^lUoo(n) or \\u - u|| L oo(n) < \\<p - <p\\L°°(<a) 

respectively. 



4 Lewy-Stampacchia inequalities and its consequences 

Let X be a real reflexive Banach space, which is a lattice with respect to a partial order 
" < ", and V is a sublattice of X (i.e. V contains the sup and inf of its elements). Then for 
every n£l one can write v = v + — v~ , where v + = v V 0, v~ = — v A are the positive and 
negative parts of v respectively. Thus, X is generated by the cone P, P — P = X, 

P = { v e X : v>0}. 

The subspace of V (the dual space of V) generated by the cone 

P> = {v 1 G V : («', v) > 0, Vv G P}, 

is called the dual order of V and denoted by V*, i.e. V* = P' — P' . We shall consider an 
operator A : X — > V' with the properties: 

a) hemi-continuous, i.e. the mapping t — > (A(u + tv),w) is continuous on [0, 1], u,v G X, 

w G V; 

b) coercive, i.e. 3wq G X such that 

•r ta i- (Au,U — Wq) 

it u — Wo G V, lim — = oo, u G X, 

\\u\\— »oo U 
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c) strictly T-monotone, i.e. 

(Au - Av, (u - v) + ) > 
for all u,v G X such that / (u - v) + G V. 

We recall that L > M for L, M G V , if L - M G P', i.e. L — M is positive on the positive 
elements of V. Then, if A is strictly T-monotone, it satisfies a weak comparison principle: if 
Au < Av, then u < v. 

In this framework we recall the following interesting fact (see [25]). 

Theorem 4.1. Let A : X — > V' be a strictly T-monotone, coercive and hemi- continuous 
operator. Let also u,v G V be such that Au,Av G V*. Then A{u A v), A(u V») £ V* and 
A(u Aw) > Au A Av, A{u V v) < Au V Av for the dual order in V' . 
Consequently, if Aui G V* , then for Vi G I C {1, 2, . . . , iV} 



In fact, we can take X = W 1 ' ^), V = W^ G (n) and the operator A defined by (3.2) 
under our assumptions. 

Now we prove the Lewy-Stampacchia inequalities in abstract form. 

Theorem 4.2. Let X be a real reflexive ordered Banach space, V be a closed subspace of 
X, which is a sublattice of X , A : X — > V be a strictly T-monotone, coercive and hemi- 
continuous operator. Let also there are given two elements <p,ip G X, -ip < <p, and L G V . If 




(a) 



3 A G V such that A > L and A > Aip in V' 
3 A G V' such that A < L and A < Aip in V 1 , 



(b) 



(tp — v) + G V and (v - (p) + G V, Vv G V, 



and let u be the solution of 



BElJ: (Au — L,v — u) > 0, Vi> G K% 



(4.1) 



where 



K^ = {ueV <u<p}. 



(4.2) 



T/ien t/ie following dual estimates hold 



\<Au<A in V'. 



(4.3) 



In particular, if L, Aip, Aip G V*, f/ien also Au G F* ; and (4.3) gives 



L A Aip < Au < LV Aip in V* . 



(4.4) 
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Proof. The assertion (4.4) assertion holds, since, if L, Aip, Aip G V* , one can take A = 
L A Aip = L — (L — Aip) + and A = L V Aip = L + (Aip - L)+ in (4.3). 

To prove the upper bound of (4.3), consider the unique solution z G V of auxiliary variational 
inequality 

z < u : (Az - A, w - z) > 0, Vw; G V, u> < u. (4.5) 

It is enough to show, that z = u, since then taking w = u — v in (4.5) for an arbitrary v > 0, 
it readily follows 

Aw - A = Az - A < in V' . 

To prove that z = u, let us first prove that z > ip. Recalling the first condition of (b) and 
taking w = z + (ip — z) + = ip V z < u in (4.5), we get 

(A - Az, (ip - z) + ) <0. 

Hence, since A > Aip, one obtains 

(Aip — A,(ip — z) + ) + (A - Az, (ip - z) + ) < 

which, by the strict T-monotonicity of A, implies, that (ip — z) + = 0. This means that z >ip. 
Let us now prove that z > u. Since u solves (4.1) and in (4.5) z < u, one has that z < (p. 
In other words, z G K^. Take w = zVn = z+(u — z) + < u in (4.5) and u = it A z = 
u — [u — z) + > ip in (4.1). Then, by addition and recalling the fact A > L from the conditions 
(a), we conclude 

(Au - Az, (u - z) + ) < (L - A, (u - z) + ) < 0, 

and, since A is strictly T- monotone, it follows that (u — z) + =0 and so z > u. 

So, the unique solution u of (4.1) is also the unique solution of (4.5). We already know, that 

this implies the upper bound in (4.3). 

To prove the lower bound in (4.3), consider the unique solution z G V of auxiliary 
variational inequality 

z > u : (Az - A, w - z) > 0, \/w G V, w > u. (4.6) 

The steps are similar as above. It is enough to show, that z = u, since then taking w = u + v 
in (4.6) for an arbitrary v > 0, it readily follows 

A - Au = A - Az < in V'. 

To prove that z = u, let us first prove that z < ip. Recalling the second condition in (6) and 
taking w = z — (z — ip) + = z A if > u in (4.6), we get 

(Az- A, (z -</?)+) < 0. 

Hence, since A < Aip, one obtains 

(A - Aip, (z - ip) + ) + (Az - A, (z - ip) + ) < 0, 

which, by the strict T-monotonicity of A, implies, that [z — ip) + =0, i.e. z < ip. 

Let us now prove that z < u. Since u solves (4.1) and in (4.6) z > u, one has that z > ip. 

In other words, z G K^. Take w = z/\u = z — (z — u) + > u in (4.6) and v = u V z = 
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u + (z — u) + < p in (4.1). Then, by subtracting and recalling that A < L from (a), we 
conclude 

(Az - Au, (z - u) + ) < (A - L, (z - u) + ) < 0, 

and, since A is strictly T-monotone, it follows, that (z — u) + = and so z < u. 

So, the unique solution u of (4.1) is also the unique solution of (4.5) (we already know, 
that this implies the upper bound in (4.3)) and the unique solution of (4.6), which, as we 
see, implies the lower bound in (4.3). ■ 

Remark 4.3. Theorem 4-2 is still true when X is just an ordered Banach space and the 
operator A is strictly T-monotone, provided the problem (4.1) ; (4.2) has solution. This result 
extends Theorem 4-1 of [25] that was restricted to the abstract one obstacle problem, and 
Theorem 4-35 of [30] was stated only for linear second order operators in Sobolev spaces, 
which was recently extended to the p-Laplacian in [28] and to general Leray-Lions operators 
in [24]. 

Let u be the unique solution of (4.1), u be the unique solution of (4.1) in K^, = {v € 
V : ip < v} (i.e. u is the solution of one (lower) obstacle problem) and let u be the unique 
solution of (4.1) in = {v € V : v < ip] (i.e. u is the solution of one (upper) obstacle 
problem). We can recover the Lewy-Stampacchia inequalities for the one obstacle problem 
easily in the following way. 

Proposition 4.4. i) If Aip > L, then u = u, and one can rewrite (4.4) as 

L < Au <L\I Aip in V* . 
ii) If Aip < L, then u = u, and one can rewrite (4.4) as 

L A Aip < Au < L in V*. 
Proof, i) By taking v = u — (« — ip) + = u A ip € Ky, in the lower-obstacle problem 
u € IK,/, : (Au — L, v — u) > 0, Vv £ K.^, 

we get 

{Au — L, {u-p) + ) < 0. 
On the other hand, by assumption, 

(L-Aip, (u- ip) + ) < 0. 

Recalling the strictly T-monotonicity of A, by addition, we conclude (u — (p) + = and so 
u<ip, which means, that u 6 K^. The uniqueness of the solution gives u = u. 

ii) Similarly, by taking v = u+ (ip — u) =u\/ ip £ in the upper-obstacle problem 

u G K v : (Au - L, v - u) > 0, VU G K^, 

we get, using 

(Aip - Au, (V> - u) + ) = (Aip — L,(tp — u) + ) + (L - Au, (V> - u) + ) < 0. 

So (ip — u) + = 0, which means that u G K^, and we conclude by the uniqueness of the 
solution that u = u.M 
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We have already observed that the nonlinear operator A given by (1.1) under the general- 
ized condition (1.5) of Ladyzhenskaya and Uraltseva type satisfies the hypothesis of Theorem 

1 G 

3.2 in Wq' (Q). In particular, if we also assume the existence of C > 0, such that 

< C, A^<C and Atp > -C a.e. i£fi, (4.7) 

the solution of (3.3) satisfies 

-C < (/ A A<p)(x) < Au(x) < (/ V Aip){x) < C a.e. x E CI, 
and we conclude that 

Au G L°°(fl). 

Hence, the regularity of the solution of the two obstacles problem (and, similarly, for each 
one obstacle problem as well) is the same as bounded solutions of the respective equation 
without obstacles, as it was established in [20]. So if, we assume, in addition to the continuous 
differentiability in t of a(x,t), that for some constants Cp > and < (3 < 1 

\a(x,t)t - a(y,t)t\ < Cp(l + a(x,t)t)\x - yf, Vx,y G H, t > 0, (4.8) 

from [19], [20] we immediately conclude the following interesting regularity result (which for 
the p(x)-Laplacian operator was re-established in [13]). 

Theorem 4.5. Under the assumptions (1.5), (3.4), (4.7) and (4.8) the solution u of the 
obstacle problem (3.3) with <p, tp G L°°(ft) is C^fi), for some < a < 1. If also dtt G C 1 '", 
i/ien « G C 1 '"^). 

Remark 4.6. The Holder continuity of the gradient to the one or two obstacles problem has 
been obtained by Lieberman [19] when the obstacles are also C l,a . However, the condition 
(4.7) on ip and ip do not imply that cp, ip belong also to C 1,a . So the result of Theorem 4-5 
does not follow nor imply Lieberman's C 1,a regularity results. 

5 Applications to Systems 

5.1 N-membranes problem 

For the operator A the N-membranes problem consists of: to find u = (ui, U2, ■■■,un) G 
satisfying 

N N 

V] / a(x, \Vui\)Vui ■ V(vi - m)dx > V] / fi{vi - Uijdx, V(ui, ...v N ) G K N , (5.1) 
i=i Jn i=i Jn 

1 G Af 

where is the convex subset of the Orlicz-Sobolev space [W^{n)] N , defined by 

K N = {( Vl ,...,v N ) G [Wo' G (n)] N : v x > ... > v N a.e. in 0}, 

and /i,...,/jv G L° (fi). As in Theorem 3.2 the existence and uniqueness of the solution to 
(5.1) follows easily. 
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Theorem 5.1. If a(x,t) satisfies to the conditions above (namely (1.6)), then the solution 
(ui, ujv) of the N -membranes problem for A satisfies the following Lewy-Stampacchia type 
estimates: 

fi A Au\ < ft V ... V/jv 
/l A / 2 < Au 2 < / 2 V ... V 

fl A ... A /jv-1 < AujV_i < /at-i V f N 
fi A ... A f N < Au N < f N 

a.e. in Q. 

Proof. Observe that choosing (v,u 2 , ...un) G IK/v, with u G K U2 , we see that ui G K U2 
solves the "lower-obstacle problem" with f = fi, and so (recall Proposition 4.4) 

/i < < /i V ^4ti2 a.e. in fl 

Analogously, we see that Uj G solves the two-obstacles problem with f = fj, j = 

2,3, ...TV - 1, and satisfies, by (4.4), 

fj A Auj-i < Auj < fj V Auj + i a.e. in Q. 

Since ujy G IfC"^- 1 , then (recall Proposition 4.4) 

fN A Aun-i < Au^y < /at a.e. in $7. 

The proof concludes by simple iteration. ■ 

Since for each i = 1, 2, . . . , N, the respective Aui maybe controlled in L°°(Q) by 

A fj<A Ui < v ./•, 

l<j<i i<j<N 

if a(x, t) is Holder continuous in x and continuously differentiable in t, by the regularity of 
[20] we also have 

Theorem 5.2. If ft G L°°(0) (t = 1, 2, . . . N) and (1.5), (4.8) feoW, iften i/ie solution u of 
the N-membranes problem has C 1,Q! (f2) regularity. If also d£l G C 1,a , i/ten £/ie solution u to 
N-membranes problem belongs to [C 1,a (Jl)\ N for some < a < 1. 

As in [5] we may also approximate the variational inequality using bounded penalization. 
Defining 

c f fl H r- /i . \ 

£ = max <^ : % = 1, . . . , N \ , 

Ci=i^-(fi + --- + ft)loTi = l,...,N, 

we observe that 

/ & > 0, if i > 1, 

I - 6-2) - & - &-i) = /< - /i-i, if * > 2. 
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For £ > 0, let 9 £ be define as follows: 



6 £ : R -> M, s i-> < 



' 0, if s > 0, 
s/e, if — e < s < 0, 
— 1, if s < —£. 



1 G 

The approximate problem is given by the system for u\ 6 W ' (fi) 

Ait? + £A« - <+i) - &-i0 £ «_i - <) = /i, m i = 1, . . • ,N, (5.2) 

with the convention Uq = +oo, u £ N+l = — oo. Since the operator A is strictly T-monotone, 
then arguing as in [5], we get: 

Theorem 5.3. If the operator A satisfies the assumption 0-*>\> then 

(i) the problem (5.2) has a unique solution (u\, . . . , u £ N ) G [W ' (^)] 7V - TTiis solution satisfies 

<<<_! + £ /or » = 2,..., AT. 

(ii) («f, . . . , -ufy) — > (tii, . . . , ujv) ^ \Wo' (^)] N strongly as £ — > 0, where (u\, . . . , ujv) is the 
solution of the N-membranes problem. Under assumptions of Theorem 5.2 the convergence 
holds in [C 1 ^)]^. 

Proof. As observed in [5], the bounded penalization B £ 

N 

{B £ v, w) = V / [ii9 £ {vi - v i+ i) - &_i0 e (ui_i - Wi)]^i 

is monotone in [L 1 ($7)] Ar . Therefore (i) follows exactly in the same way as in Proposition 2.1 
of [5]. 

Since (B £ \, v — u e ) = for any v £ K/v, the weak convergence of (ii) follows also exactly as 
in [5] by monotonicity arguments, since we have 

(Au £ ,v-u £ ) > (L,v-u e ), VvgKat. (5.3) 

Here we have used the notations 

N N , 

(Aw,v) = V(Awi,t>i) and (L,v) = V / fm dx. (5.4) 

1=1. 1=1 J ^ 

Since u £ — ^ u € Kjv weakly in [W^ (£l)] N , setting (5.2) in variational form, using (5.3) with 
v = u and taking the limsup we obtain the strong convergence (using Lemma 3.5). Finally 
if also f € {L°°(£1)} N then the penalization term in (5.2) is also bounded in L°° and the u\ 
are uniformly bounded in C l ' a M 

5.2 A Quasi- Variational Problem 

Another interesting application of the Lewy-Stampacchia inequalities is its application in 
studying quasi-variational inequalities. Some of these problems are related to a stochastic 
switching game. 

In this section we assume that dQ £ C 1,13 for some < /3 < 1, and the operator A has 
the C l ' a regularity property. 
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We consider the following quasi-variational problem (with the notations (5.4)): 
find u G K(u), such that, 



(Au - L, v - u) > 0, Vv € K(u), (5.5) 



where 



K(u) = {v G [W 1,G (fi)] N : *i(u) < Wi < ^(u), i = 1, 2, . . . , N}, 
f = (/i,... ,/jv) G [L g (0)] w , and for v = . . . ,ujv) we set 

^i(v) := /\( w j + Vij) (Vij are positive constants), 

VPj(v) := \J (vj — ipij) (tpij are positive constants), 

for i,j = l,...,N, i.e. we consider the problem (5.1) with Kjy = K(u) (the obstacles them- 
selves depend on the solution). A similar problem for linear operators (and for the one 
obstacle) was considered in [14] and [30]. 
We set 

N N 

a* : = A fc> and v := V f* ( 5 - 6 ) 

i=i i=i 

and we consider the unique solution u° and if in W 1 ' (fi) solving, respectively, 

Au° = ijl and An = v in (5.7) 

Assuming G [L 00 ^)]^, we have «°, «° G C l,a (£l) (see [20]) and we may define the constant 
Ao by 

Ao = maxH — minn > 0. 
n n 

Theorem 5.4. Let ft G [L°°(n)] N and suppose 

fij + Ak > A > 0, Vi,j,k = 1, . . . , iV. 
T/ten (5.5) admits at least a maximal solution u and a minimal solution u. 
Proof. Let /x, v G L°°(f2) be given by (5.6). By regularity, the set 

B := {v = G [Wq 1 ' ^)]^ fi < Avi < v in Q, i = l,...,N} 

is bounded in [C 1 ' a (r2)] iV for some a > 0. By comparison we have vP < V{ < vP in Q, 
i = 1,...,N, for every v G B and, of course, u° = (u°, . . . , u°) and u° = (vP, . . . , u° also 
belong to B. Therefore, if v G B we have 

maxfj — mintifc < Ao, Vj, k 
n a 

and 

K(v) = { z = ( Zl ,...,z N )e [W^ G (Q)] N : *i(v) < z % < $i(v), i = 1, . . . , iV} ^ 0. 
If we denote by w = cr(v) the unique solution of 

wGK(v): (Aw-L,z- w) > 0, Vz G K(v), (5.8) 
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by Lewy-Stampacchia inequalities and Theorem 4.1, it satisfies 

( f\ AvA Afi<A Wi <(\/ AvA V fi in a 

Consequently, w G B and cj(B) C B. 

Using (5.8) we define u 1 = cr(u°) and by the weak maximum principle, we observe that 
u\ < v° in f2, for all i = 1,. . . ,N. By iteration, let u m+1 = a(u m ) for m = 1, 2, . . . and 
observe that, if componentwise uf 1 < * nen 

^(u m ) < ^(u"^ 1 ) and $,(if) < $i(u m4 ). 

By monotonicity with respect to the obstacles, we conclude then 

< _™ j = i 5 . . . 5 # for all m = 0, 1, 2, . . . 

Therefore we have constructed a decreasing (and bounded in [C 1 ' a (fi)] iV ') sequence u™ in D, 
which converges as m — > oo to u G B in [C 1 (0)] 7V . In particular, this implies, u 6 K(u) and 

$(11™) -> $(u) and *(u m ) -> *(u) in [Vy 1 ' G (0)] Jv . 

Since u m+1 G K(u rn ) also 

^_ m+ l L _ _ m+ i^ > Q) v _ G K( ^my 

By Theorem 3.6, u solves (5.5). 

Analogously, for m = 0, 1, 2, . . . we construct an increasing sequence 

u m < u m+1 = <r(u m ) -> u in [C 1 ^)]^ n B 

and easily conclude that u G K(u) also solves (5.5). 

If u = (ni, . . . , un) is any other solution to (5.5), by the Lewy-Stampacchia inequalities 
we have for alH = 1, . . . , N 




and so, by comparison, u < m < Tr. Hence u] < Ui < u\, by monotonicity of and <I>, and 
so also u™ <Ui< u^ 1 for all m by recurrence. 

We then conclude u { < Ui < u-i, i = l,...,iV that yields the minimality of u and 
maximality of u.M 

Remark 5.5. The uniqueness of the solution to (5.5) is an open problem. 

References 

[1] Adams R.A, Sobolev Spaces, Academic Press [A subsidiary of Harcourt Brace Jovanovich, 
Publishers], New York-London, 1975, Pure and Applied Mathematics, Vol. 65, xviii+268 
pages. 

[2] Acerbi E., Mingione G., Regularity results for a class of functionals with non-standard 
growth, Arch. Ration. Mech. Anal. 156 (2001), no. 2, 121-140. 



16 



[3] Antontsev S.N., Rodrigues J.F., On stationary thermo-rheological viscous flows, Ann. 
Univ. Ferrara, Sez., Sci. Mat. 52 (2006) pp. 19-36. 

[4] Antontsev S.N. and Shmarev S.I., Elliptic equations with anisotropic nonlinearity and 
nonstandard growth conditions, Handbook of Differential Equations: stationary partial 
differential equations (M. Chipot and P. Quittner, eds.), vol. 3, Elsevier (2006) pp. 1- 
100. 

[5] Azevedo A., Rodrigues J.F., Santos L., The N-membranes problem for quasi-linear de- 
generate systems, Interfaces and Free Boundaries 7 (2005), pp. 319-337. 

[6] Benkirane A., Elmahi A., An existence theorem for a strong nonlinear elliptic problem in 
Orlicz spaces, Nonlinear Analysis, No. 36 (1999), pp. 11-24. 

[7] Challal S., Lyaghfouri A., Holder Continuity of Solutions to the A-Laplace Equation 
Involving Measures, Communications on Pure and Applied Analysis, vol. 8, No 5, (2009), 
pp. 1577-1583. 

[8] Chen Y., Levine S., Rao M. Variable exponent, linear growth functionals in image restora- 
tion, SIAM J. Appl. Math. 66 (2006), no. 4, pp. 1383-1406. 

[9] Diening L., Theorical and numerical results for electrorheological fluids, PhD Thesis, 
University of Freiburg, 2002, Ph. D. Thesis. 

[10] Diening L., Maximal function on Musielak- Orlicz spaces and generalized Lebesgue spaces, 
Bull. Sci. Math. 129 (2005), no. 8, p. 657-700. 

[11] Donaldson T.K., Trudinger N.S., Orlicz- Sobolev Spaces and Impedding Theorems, Jour- 
nal of Functional Analysis 8 (1971), 52-75. 

[12] Elmahi A., Meskine D., Elliptic inequalities with lower order terms and L 1 data in Orlicz 
spaces, J. Math. Anal. Appl. 328 (2007), pp. 1417-1434. 

[13] Fan X. Global C 1,a regularity for variable exponent elliptic equations in divergence form, 
J. Differential Equations 235, no. 2, 397-417 (2007). 

[14] Joly J., Mosco U. A propos de I 'existence et de la regularite des solutions de certaines 
inequations quasi-variationelles, J. Functional Anal. 34 (1979), pp. 107-137. 

[15] Gossez J., Nonlinear elliptic boundary value problems for equations with rapidly (or 
slowly) increasing coefficients, Trans. Amer. Math. Soc. 190 (1974), 163-205. 

[16] Gossez J., A strongly nonlinear elliptic problems in Orlicz- Sobolev spaces, Proc Sympos, 
Pure Math., 45. Part 1, (1986), pp. 455-461. 

[17] Gossez J., Mustonen V., Variational inequalities in Orlicz- Sobolev spaces, Nonlinear 
Anal. Theory Meth. Appl. , 11 (1987) pp. 379-392. 

[18] Kufner A., John O., Fucik S., Function spaces, Noordhorff, Leyden, 1997. 

[19] Lieberman G.M., Regularity of Solutions to Some Degenerate Double Obstacle Problems, 
Indiana University Mathematics Journal, Vol. 40, No. 3 (1991), pp. 1009-1028. 



17 



[20] Lieberman G.M., The Natural Generalization of the Natural Conditions of Ladyzhen- 
skaya and Uraltseva for Elliptic Equations, Commun. in Partial Differential Equations, 
16 (2 and 3) (1991), pp. 311-361. 

[21] Le V.K., Schmitt K., Quasilinear elliptic equations and inequalities with rapidly growing 
coefficients, J. London Math. Soc. (2) 62 (2000), pp. 852-872. 

[22] Martinez S., Wolanski N., A minimum problem with free boundary in Orlicz spaces, 
Advances in Mathematics 218 (2008), pp. 1914-1971. 

[23] Mihdilescu M., Radulescu V., Neumann problems associated to nonhomogeneous differ- 
ential operators in Orlicz- Sobolev spaces, Annalies de l'institut Fourier, tome 58, no. 6 
(2008), p. 2087-2111. 

[24] Mokrane A., Murat F., The Lewy-Stampacchia inequality for bilateral problems, Ricerche 
Mat. 53 (2004) pp. 139-182. 

[25] Mosco U., Implicit Variational Problems and Quasi Variational Inequalities, Lecture 
notes in Mathematics, No 543, Springer Verlag, Bruxells, 1976, pp. 83-156. 

[26] Musielak J., Orlicz Spaces and modular spaces, Lecture Notes in Mathematics, vol. 1034, 
Springer- Verlag, Berlin, 1983, iii+222 pages. 

[27] Rodrigues J.F., Obstacle Problems in Mathematical Physics, North-Holland Mathemat- 
ics Studies 134 (Notas de Matematica 114), Elsevier Science Publishers B.V., 1991. 

[28] Rodrigues J.F., Stability remarks to the obstacle problem for p-Laplacian type equations, 
Calc. Var. 23 (2005), pp. 51-65. 

[29] Ruzicka M., Electrorheological Fluids: Modeling and Mathematical Theory, Lecture 
Notes in Mathematics (11748, Springer, Berlin, 2000. 

[30] Trioaniello G.M., Elliptic Differential Equations and Obstacle Problems, 1987, Plenum 
Press, New York. 

[31] Yamada N., Viscosity Solutions for a System of Inequalities with Bilateral Obstacles, 
Funkcial. Ekvac, 30 (1987), pp. 417-425. 

[32] Zhikov V.V., Averaging of functionals of the calculus of variations and elasticity theory, 
Izv. Akad. Nauk SSSR, Ser. Mat. 50 (1986), no. 4, p. 675-710,877. 

[33] Zhikov V.V., Meyer-type estimates for solving the nonlinear Stokes system, Diff. Uravn. 
33 (1997), no. 1, p. 107-114, 143. 

[34] Zhikov V.V., On the Technique for Passing to the Limit in Nonlinear Elliptic Equations, 
Functional Analysis and Its Applications, Vol. 43, No. 2, pp. 96-112, 2009. 



18 



